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ABSTRACT 



An analysis of jet stream structure and energetics is performed 
utilizing an objective scheme that calculates wind components normal 
and parallel to the jet axis, with provisions to explicitly calculate 
the eddy components of variables. Most of the meteorological data used 
is generated by a diagnostic-balanced, e-equation model. Results indi- 
cate that above 500 mb observed winds and temperatures should be used 
versus computer generated data. 

There is seme indication that between the jet entrance and exit 
regions the mean transverse vertical circulation reverses from thermally 
direct to thermally indirect. Kinetic energy is imported from below 
the jet into the jet stream region primarily by the mean pressure 
interaction at 500 mb. 

Thermal balance calculations indicate that the mean horizontal 
temperature advection and vertical transverse motion are primarily 
responsible for maintenance of the temperature structure about a jet. 
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area enclosed by lateral boundaries 
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performed along lateral boundaries 

velocity component normal to the lateral 
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I. INTRODUCTION 



Since the mid-1940's, when the existence of jet streams in the 
troposphere was documented by Rossby and Staff (1947) , many studies of 
the Polar Front Jet Stream have been conducted. Notably, these studies 
have yielded a considerable amount of empirical knowledge about the 
Polar Front Jet. Unfortunately, the energetics of this jet are less 
well understood. 

One important reason for this lack of understanding has been 
pointed out by Mahlman (1970) . Considerable difficulty exists in 
defining an adequate curvilinear coordinate system for such a strongly 
time-dependent phenomenon . Seme investigators have attempted to over- 
come this difficulty and have conducted studies of jet streams in 
curvilinear coordinates, with most interesting results . Riehl and 
Fultz (1957, 1958) in their steady, three-wave rotating dishpan experi- 
ment, investigated the dishpan jet stream in curvilinear coordinates, 
centered on the jet axis. They found net ascent poleward and net 
descent equatorward of the strong westerlies when averaging was per- 
formed around latitude circles. However, when the averaging was per- 
formed along the jet axis, the vertical circulation was reversed. 

This implied that the mean transverse circulation about the dishpan 
jet was thermally direct, while the mean meridional circulation in 
the same latitude operated in a thermally indirect sense. This 
suggested that the mechanisms responsible for maintaining the jet 
streams in a geographic coordinate system may be different from those 
in a curvilinear coordinate system. 



9 



The most detailed study of an atmospheric jet stream in curvilinear 
coordinates was performed by Krishnamurti (1961) , on the subtropical 
jet stream, covering the period Decanter 1955 and January-February 
1956. The subtropical jet stream is particularly suitable to analysis 
in curvilinear coordinates since its day-to-day variability is rather 
small, especially when compared with the polar front jet. Krishna- 
murti's study also showed the existence of a thermally direct trans- 
verse circulation about the jet core. This agreed with the net 
transverse circulation about the jet in curvilinear coordinates. . The 
mean meridional circulation in the latitudes of the subtropical jet 
operates in the same sense as the transverse circulation, in contrast 
to the dishpan jet. 

The preceding studies avoided the time-dependent problem associated 
with jet streams. In the case of the dishpan jet the flew was in 
steady state with only longitudinal translation, and in the subtropical 
jet case the flow patterns were nearly stationary. In addition, both 
studies were limited to sane extent, in that no explicit calculations 
of eddies in a curvilinear coordinate system were made. 

Until recently no attempt has been made to extensively investigate 
the Polar Front Jet Stream in curvilinear coordinates since it is the 
most time-dependent of the major tropospheric jets. Mahlman (1970) 
minimized this time-dependent problem by choosing a Polar Front Jet 
where the time variations were not particularly large. While his 
study was quite extensive, it was limited primarily by two factors. 
First, an excessive amount of time was required to define the grid 
systan, extract data, and perform the necessary calculations in the 
curvilinear coordinates. The second difficulty was in explicitly 
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computing eddy components of the wind in a curvilinear coordinate system. 
In this system an eddy quantity was defined as the value of a variable 
at a given point minus its mean value along the same line parallel to 
the jet axis. As Mahlman (1970) pointed out, in a curvilinear coordi- 
nate system, the system tends to be aligned nearly parallel to the flew. 
Thus, the conputation of wind eddies was extremely sensitive to the 
orientation of the jet axis and small variations in the wind flow. 

The first intent of this study was to develop a method for rapid, 
objective analysis of a jet stream. Meteorological input data was 
generated by a diagnostic-balanced, co-equation program (Krishnamurti 
1968), which required only geopotential heights. Thus both the 
divergent and nondivergent components of the wind are included in the 
study. Some discussion of the advantages and disadvantages of the data 
used, plus seme steps necessary to minimize the disadvantages, is 
contained in Section II-B. 

The objective computational scheme defined a coordinate along a 
jet stream axis and coordinate lines at distances 2.5 and 5.0 degrees 
latitude normal to the jet axis on either side (Fig. 1) . Then the 
shceme linearly interpolated data from the various fields; co,T,$,U and 
V components of the total balanced winds. The scheme calculated wind 
components normal and parallel to the jet axis. Line or area averaging, 
where appropriate, was performed to calculate the eddy components of 
variables. A detailed description of this scheme is contained in 
Section II-D. 

Once the above computational scheme was developed, several checks 
were performed. First, results of the data interpolated from the 
initial data fields by the computational scheme were compared with 
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independent interpolations done by hand. Second, using the same time 
periods and jet axes at 300 mb that Mahlman (1970) used, the jet stream 
kinetic energy balance was calculated. The basic data fields and com- 
putational procedures were discussed by Mahlman (1970) . The results 
obtained were then compared with Mahlman' s results based on subjective 
interpolations and hand calculations. 

The results of these checks were quite successful. Comparison of 
the data interpolated from the initial data fields by the computer and 
by hand showed very small differences, especially when averaged over 
the five time periods involved. This meant that the data interpolated 
by the objective scheme could be used with confidence in calculating 
the jet stream kinetic energy balance. Comparison of Mahlman' s and the 
objective scheme calculations for the jet stream kinetic energy 
balance also shewed good results which will be presented in Section IV. 

Having satisfactorily developed the objective computational scheme, 
three further objectives for this research were set. First, due to 
certain inadeguecies in the available data, which will be discussed in 
detail in Section II-B, jet axes at 200 mb were defined. The jet 
stream kinetic energy balance with respect to these axes was then cal- 
culated as it was for the axes at 300 mb mentioned previously. Second, 
since Mahlman' s results were based on five time-period averages, it 
was felt that it would be of seme interest to study a single time 
period. Therefore, the jet stream for 16 November 1200 OVir, 1966, 
was analyzed in detail. This period appeared to present the most 
explicit situation of a jet stream maximum of the five time periods 
studied. 
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Lastly, Mahlman (1970) discussed the role of eddy heat flux in the 
maintenance of the temperature structure about the jet. Mahlman's 
analysis indicated that eddy heat flux might counteract the expansion 
cooling and compression heating caused by the mean circulation about 
the jet. Therefore, the third objective of this research was to 
explicitly investigate the role of eddy heat flux in the maintenance 
of the temperature structure about a jet stream. 

II. RESEARCH PROCEDURES 

A. DESCRIPTION OF THE JET STREAM COORDINATE SYSTEM 

The jet stream coordinate system for this project was defined, 
consistent with Mahlman (1970) , in the following manner. A central 
coordinate was established, in this case, along the jet stream axis. 
Parallel coordinate lines were defined at distances 2.5 and 5.0 degrees 
latitude, normal to the jet axis in either direction. In an (s, n, p) 
curvilinear coordinate system, s and n were respectively the coordinates 
parallel and normal to the jet axis, and £ was taken as the vertical 
coordinate. Separation between grid points along the £ coordinates, 
other than the jet axis, was not uniform. This was the major difference 
between this grid and the one Mahlman used with uniform separation along 
all £ coordinates. Once this jet stream-oriented coordinate systan was 
defined by specifying an axis at a single pressure level, it is consistent 
through all levels in the vertical (See Fig. 1) . 

B. MEIEOROIOGICAL DATA USED 

The case selected for analysis was an area centered over continental 
United States as described by Mahlman (1970) . Data at five synoptic 
observation times was used, frcm 15 November 1200 QMT to 17 November 
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Fig:, l Curvilinear grid for jet stream of 15 November I9b6, 
1200 GMT 




Fig. 2 Geopotential heishts(m) at 300 mb on 15 November 
1966, 1200 GMT 
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1200 GMT, 1966. During this period the flow at jet stream level was 
anticyclonic between troughs off the west coast and over Hudson Bay in 
the east. Although there was scare change in the location of the jet 
stream, no significant development or large-scale syston movanent 
occurred. The jet stream axes did shew a tendency to shift to the 
southeast with the jet maxima traveling eastward along the axis. Thus, 
this case appeared to present a satisfactory steadiness of the jet 
axis for analysis (See Fig. 2) . 

As mentioned in the introduction, all' meteorological data in this 
research were generated by the diagnostic balance, co-equation model, 
developed originally by Krishnamurti (1966, 1968) . In this model, the 
equations are scaled for motions of characteristic Rossby number less 
than unit (R q < 1) . Initial input data was obtained fran the National 
Meteorological Center's objectively analyzed geopotential fields at 
1000, 850, 700, 500, 300, and 200 mb. The stream function (i|>) was 
related to the geopotential field through the non-linear balance 
equation, 



V£V(rt-2J(|i,||)=v2* 



(1) 



If the ellipticity condition for Eqn. (1) 



V 2 <H-l/2f 2 -Vf*V4i>0 



( 2 ) 



is not satisfied, then Eqn. (1) is replaced by 



V-fVi|p=V 2 $-2J (u fi ,V 6 ) 



(3) 
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The meteorological data output was in staggered levels as shewn in 
Fig. (5). Fields of ip, $, and appear at the 1000, 800, 600, 400, 
and 200 mb levels, while up T, and 0 are given at the 900, 700, 500, 

300 mb levels. The horizontal grid points were spaced 2.5 degrees 
latitude and longitude apart. The grid extends frem 70W to 135W and 
frem 25N to 60N. Six additional grid points in the zonal direction did 
not contain initial data, but were used to provide cyclic continuity 
for any given dependent variable. These grid point values were inter- 
polated frem values at the first two and ' last two of the grid points 
with real initial data. There was no staggering of variables in the 
horizontal. For a concise description of the computational procedures 
of this model, see Krishnamurti (1968) and Mahiman (1970) . 

Seme limitations of the data have appeared. First, the jet stream 
chosen, while of moderate intensity, was not particularly well-defined. 
For most of the time periods, two distinct jet axes, whose separation 
was less than 10 degrees latitude, were present. During the last two 
time periods, 17 November 0000 GMT and 17 November 1200 GMT, 1966, a 
moderate jet stream associated with the lew center over Hudson Bay 
impinged on the northeast comer of the grid, merging with the jet 
stream being analyzed. This caused sane difficulty in defining axes 
and in interpreting the results. The manner in which this problem was 
handled is discussed in Section II-C. 

Second, the ellipticity problem of the balance equation along with 
erroneous geopotential heights adjustments had to be considered. When 
the ellipticity condition for Eqn. (1) was not satisfied, Eqn. (3) was 
substituted for Eqn. (1) . This adjustment was frequently required on 
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the anticyclonic side of the jet stream. Ellsaesser (1968) has can- 
pared the calculated winds of the balance equation with the observed 
winds. In places he found the mean jet stream at 300 mb was seventy 
percent stronger than indicated by the balance equation. Calculated 
wind speeds just equatorward of the jet stream were also found to be 
significantly lower. Calculations bear out this disadvantage of the 
non-1 inear i zed balance equation. 

Errors in data handling also appeared. This problon was not 
particularly serious, since it was obvious when the error was large. 
Subjective adjustments to the geopotential heights were made when 
necessary. 

Since it was not an objective of this research to analyze the 
balance equation irodel, no attempt was made to change the model to 
minimize the ellipticity problem. However, when interpreting the 
results this limitation was kept in mind. 

Third, having the vertical motion defined at different levels than 
the kinetic energy and temperature is an obvious limitation. The 
vertical advection terms were calculated by linearly interpolating the 
values between levels. While not considered serious this prcblan must 
be kept in mind when interpreting the results. 

C. DESCRIPTION OF JET STREAM AXES ANALYZED 

For descriptive purposes the jet axes analyzed were divided into 
two groups, those based on observed winds and those based on the com- 
puted (balanced) winds. The first group used jet axes defined frcm the 
National Meteorological Center's maps of observed 300 mb winds and 
those were the ones used to compare with Mahlman's results. 
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Fig. 3 



Jet stream axes at 300 mb, defined from 
the observed winds 




Fig. 4 



Jet stream axes at 200 mb, defined from 
the total balanced winds 
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The second group used jet axes at 200 mb defined frcm balanced 
winds (See Fig. 4). The jet streams thus specified were of moderate 
strength, fairly well-defined, and presented a reasonable representation 
of the observed winds at 200 mb. Unfortunately, the definition of the 
jet stream axes was not particularly good east of 95W. Consequently, 
only that portion of the jet axes west of 95W was analyzed. These five 
axes gave a good picture of jet stream flow frcm SW to NE but only a 
half wave-length was being analyzed. 

D. DESCRIPTION OF OBJECTIVE COMPUTATIONAL SCHEME 

A computer program was developed for the IH4 360/67 system of the 
W. R. Church Computer Center which performed the calculations described 
in the following section. The remainder of this section is a general 
description of the computational scheme used. Primary emphasis is 
placed on operations performed and principles applied. The generality 
and flexibility of this program is also discussed in seme detail. For 
a detailed discussion of the mechanics concerning how the operations 
were performed, see Appendix A. 

The computer computational scheme was divided into four subroutines: 

1 . data input 

2. development of the curvilinear grid to be superimposed on 
the geographic data fields at all levels 

3. interpolation of data from geographic to .curvilinear system 

4. computations performed using a curvilinear coordinate system. 
The data input subroutine was simply the reading of the initial data 
fields and the definition of the jet stream axis in geographic 
coordinates . 
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The curvilinear coordinate system developed in the next subroutine 
was a simple design (See Fig. 1) . The points along the jet stream axis 
were located at two degree latitude intervals along the jet axis. Then, 
utilizing the slope of the jet axis at each point in a schene described 
in Appendix A, grid points were established at 2.5 and 5.0 degrees 
latitude normal to the axis. The grid was constant for all levels in 
the vertical. The area enclosed by the curvilinear grid and the distance 
between grid points was also calculated. 

Data was interpolated fron the geographic grid to the curvilinear 
grid in the third subroutine. A linear interpolation in two-dimensions 
was used for this transformation. After the data was interpolated fron 
geographic to curvilinear coordinate points, a linear interpolation in 
one-dimension (vertical) was performed to obtain values of co and 0 

between the staggered levels. 

In the fourth subroutine with fields of temperature , vertical 
motion, geopotential heights, and the u and v components of the total 
balanced wind at the curvilinear grid points available, the first cal- 
culation was to convert the u and v wind components to components 
parallel and normal to the jet axis at each point. This is done by the 
scheme described in Appendix A, utilizing the jet axis slope. Next, the 
eddy components of each variable were calculated. Either a line or area 
averaging scheme, as appropriate, was used to obtain the mean of a 
variable along the coordinate. Then this mean was subtracted from the 
value of the variable at each point. Both eddy and mean values were 
required for the kinetic energy balance and the eddy heat flux calcu- 
lations. To facilitate the discussion, the equations used and the 
results are discussed in Sections III and IV. 
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One of the advantages of this program was its generality and flexi- 
bility. This program was made general so that it might be used with 
other than polar-stereographic maps and 27 by 15 grid arrays. The pro- 
gram was flexible in that there is no restriction that the n coordinates 
be at 2.5 and 5.0 degrees latitude fran that jet axis, or that there be 
only two coordinate lines on either side of the jet axis. The £ 
coordinate lines were not limited to 2 degrees latitude distance along 
the axis. There may be any number of £ coordinates, and of any interval 
length desired. 



III. EQUATIONS 

In this report, the averaging notation was that used by Mahlman (1970) . 
The mean of an arbitrary quantity, Q, was given by [Q] ^ . The sub- 
script enclosed by small parentheses indicates the coordinates over 
which the averaging was performed. The deviation of Q from that average 
was given by (Q) , i.e., Q=[Q] ^ + (Q) ^ for averaging processes along 
the jet. When averaging normal to the jet, n was substituted for s. 

Consistent with Mahlman (1970) the following equation was used to 
obtain the kinetic energy balance of the jet stream. 



a_ 

at 



(s) + ^ V J-* (s) , . dp _ 3K M 

2 in; g “ at 
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-(a 1 )- 



-ia 2 >- 
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The separate terms, labeled (a) tlirough (j) have the following 
interpretation : 

(a) flux of "mean" kinetic energy through the side boundaries; 

(b) , (c) flux of "mean" kinetic energy through the lover and 

upper boundaries, respectively; 

(d) conversion of "eddy" kinetic energy to "mean" kinetic energy; 

(e) an additional term in the conversion frcm "eddy" to "mean" 
kinetic energy which arises because the Coriolis parameter (f ) 
varies along the axis over which the initial average is 
performed. It is interpreted as a conversion term because it 
appears with an opposite algebraic sign in the "eddy" kinetic 
energy equation; 

(f) mean pressure interaction term at the side boundaries arising 
frcm energy conversion in an open system; 

(g) , (h) mean pressure interaction terms at the lower and upper 

boundaries, respectively; 
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(i) conversion of potential energy into "mean" kinetic energy due 
to mean transverse circulation about the jet core; 

(j) dissipation of "mean" kinetic energy due to friction and sub- 
grid-scale mixing. 

Mahlman specified that the "mean" kinetic energy balance in the symbolism 
implied in Eqn. (4) was: 

= a, +b, + c, + f + g + h+ i - "dissipation: 

"Dissipation" was defined to be a combination of residual terms 
a 2 , b 2 , c 2 , and j. The subscripts 1 and 2 referred to the mean and 
eddy values of a term respectively. 

To investigate the role of eddy heat flux, the first law of thermo- 
dynamics was averaged along the jet axis and separated into mean and 
eddy components . 
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When discussing the results of calculations using the above 
equations, the alphabetic or numeric designations of terms are used. 
This is to simplify tables presented and to not have to continually 
repeat the meaning of terms. 
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IV. RESULTS 



To facilitate discussion, the results of this research are presented 
in the following manner. First, section IV-A presents the results of 
the ccmparision between data interpolated from the initial data field 
by the computational schane versus independent interpolation by hand. 
Second, sections IV-B and IV-C present the results of the calculation 
of the jet stream kinetic energy balance for the previously described 
axes at 300 mb and 200 mb (see Section II-C) . Section IV-B also contains 
Mahlman's (1970) energy calculation for the 300 mb jet axes 
with the results of this research. Third, in section IV-D the analysis 
of the jet stream defined at 200 mb by the total balanced winds for 16 
November 1200 GMT is presented in seme detail. Both jet stream struc- 
ture and calculations of kinetic energy are presented. Lastly, in 
section IV-E results are presented of the investigation to explicitly 
determine the role of eddy heat flux in the maintenance of the temperature 
structure about a jet stream. The conclusions drawn from these results 
are presented in section V. 

A. COMPARISON OF INTERPOIATION METHODS 

The results on the interpolation of the initial data fields from 
geographic to curvilinear grid points by both the objective and manual 
techniques are shewn in Figs. 5, 6, and 7. Fig. 5 shews the mean 
vertical motion relative to the jet axis averaged over the five time 
periods. The number above each point was obtained by the objective 
scheme, the value in parentheses was obtained by hand interpolation. 

No difference was greater than 10 x 10“ 5 mb sec - 1 , and at more than 
50 percent of the locations the difference was less than 5 x 10~ 5 mb 
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Comparison of mean vertical motion ( Lee , 10 
mb sec~' ) interpolated by the objective ‘scheme 
versus hand interpolation 
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Fig.. 6 Mean wind (knots) and temperature (deg K) 

calculated by the objective computational scheme 
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Fig. 7 Mean wind (knots) and temperature (deg K) along the 
polar front jet stream averaged from 15 November, 
1200 GMT to 17 November 1200 GMT, 1966. 

KahlmanO 970) 



i 
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sec -1 . In the mean, the vertical motion field did not agree with Riehl 
and Fultz's (1957, 1958) finding of ascent equatorward and descent pole- 
ward of the jet axis (Section I) . The first and last time periods did 
show ascent equatorward and descent poleward, but these results were 
cancelled out by ascending motion poleward of the jet axis during the 
other three time periods. 

This difference in mean values is explained by the orientation of 
the jet axes used. Fig. 3 shews that the axes used did not cover an 
entire wavelength. They were predominantly in the region of a trough 
and of a ridge. This corresponds to regions of southerly flow and 
anticyclonic curvature. Riehl and Fultz (1957, 1958) found that for 
these specific regions, ascent occupied the area. Since these two 
regions were predominant in the averaging performed, ascent should be 
found throughout, being stronger equatorward. The ascending motion was 
somewhat stronger than Mahlman (1970) showed, but the patterns were 
the same. In both cases, ascending motion through the jet core was 
found. 

Fig. 6 and 7 shew the results for the mean wind (U_) and temperature 
computed by the objective s chare and by Mahlman (1970) respectively. 

The mean winds were similar both in strength and pattern, and the mean 
temperatures were almost identical. Thus the results of this section 
indicate that the objective scheme gives comparable results to the manual 
scheme for interpolating data from geographic to curvilinear coordinate 
points. 

B. KINETIC ENERGY CALCULATIONS FOR 300 MB JET AXES 

The results for the 300 mb jet stream axes kinetic energy calcu- 
lations for both the objective scheme and Mahlman (1970) are shewn on 
Table 1. The calculations were performed from 100-900 mb over the area 
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defined by the curvilinear grid described in Section II. To acquire 
further insight, and enable more complete comparison with Mahlman's 
results, the volume was divided into regions below (300-900 mb) and 
above (100-300 mb) the jet. In addition, calculations of terms 
a^, b^, c^r d and e of Eqn. (4) are included at the bottom of the 
table. This enabled some estimate to be made of eddy terms which 
Mahlman included as "dissipation". Although the results with the 
objective scheme were generally comparable with the hand-calculated 
values, there were some significant differences. 

This was expected since the objective scheme used balanced winds 

and Mahlman used observed winds for kinetic energy computations. Where 

no horizontal wind calculations were involved, terms g and h, Eqn. (4), 

the results were almost identical. However, a definite difference 

existed, especially at upper levels, between terms where horizontal 

wind calculations were made, that is, terms a^, b^, c^, anf f in 

Eqn. (4). Thus, the problems conveming balanced winds discussed in 

section II-B must be considered vhen comparing these terms. The local 

change term ( ) best illustrated this problem, since it involved 

only u T and v T . Table 1 shows that for all layers the objective schane 
J J 

results for were lower than those calculated using observed winds . 

For this case this was in agreement with Ellsaesser's (1968) finding that 
balanced winds generally were lower than observed winds. 

The role of the flux of mean kinetic energy and pressure interaction 
were the most significant differences between the two methods of cal- 
culations. In the 100-300 mb layer Mahlman's results shew a net gain 
of +5.7 ergs atf 2 mb“ ^ec ** versus the objective schemes net loss of 
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-2.9 ergs an -2 mb~ ^sec - -'for the flux of kinetic energy, terms a^,b^, and 
in Eqn. (4). The pressure interaction terms (f, g, and h) in 
Eqn (4), both show a net gain, +6.3 versus +13.4 ergs cra -2 mb - 1 sec“ 1 . 

For the 300-900 mb layer, however, the net results are very close. 
The flux terms a^,b^, and c^, Eqn. (4) show a net loss in both cases, 
-2.9 versus -3.0 ergs an -2 mb“ 1 sec“ The pressure interaction terms 
also showed almost identical net losses, -4.4 versus -4.2 ergs 
an -2 rnb - ^ec - *. 

The large difference in the upper layer 100-300 mb, had a biasing 
effect on the total layer (100-900 mb) . For this layer the flux of 
kinetic energy was of the same sign, that is, a net loss, but the 
objective scheme's results were larger, -0.8 versus -3.0 ergs 
cm~ 2 mb” ^ec - 1 . For the pressure interaction effects, the objective 
scheme shows a net gain, +0.3 ergs an mb sec versus a net loss 
^hewn by Mahlman, -1.6 ergs an -2 mb“ ^sec - Since terms b^, c^, g and 
h, Eqn. (4), gave essentially the same results, the terms a, and f, 

Eqn. (4), in the 100-300 mb layer were the primary difference. 

These discrepancies were felt to be the result of the following 
differences. First, the objective schane used balanced versus the 
observed winds used by Mahlman. Second, the objective scheme explicitly 
calculated Vj. Thus the calculations performed were affected by the 
jet axes to the south which effected this region (Section II-B) . This 
effect was not considered by Mahlman (1970) . 

The results of the objective scheme and Mahlman' s are comparable, 
if the problems mentioned above are considered. However, it does appear 
that the results of calculations using balanced winds at or near jet 
stream level are questionable. 
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The lower portion of Table 1 shows the calculation of eddy terms 

b 2 , c 2 , d, e and thus the implied value of dissipation. Calculation 
of these terms primarily involved using u T and v T , with the previously 
mentioned problems involved. Taking these problems into account, two 
findings concerning net results appear to be significant. First, the 
net value of the eddy terms is large enough to imply that these terms, 
a 2 , b 2 , c 2 , d and e in Eqn. (4) , do have an effect on the kinetic 
energy balance. Second, the residual term j in Eqn. (4) is of sufficient 
size to suspect that dissipation of mean kinetic energy due to friction 
and sub-grid-scale mixing does occur. 

C. KINETIC ENERGY CALCUIATIONS FOR 200 MB JET AXES 

The results of the kinetic energy calculations using jet axes at 
200 mb are shown in Table 2. As in the previous section, results are 
presented for a total layer and two sub-layers. The sub-layers, however, 
are different from those shewn in Table 1. There were two basic reasons 
for this. First, the jet axes were at different levels, 200 mb versus 
300 mb. Since the balance wind data was available only to 200 mb, it 
was impossible to calculate the wind in the region above 200 mb. Second, 
from the results of the comparison in section IV-B and the problems 
involved with the data (Section II-B) , it was felt a different division 
of sub-layers was called for. Thus, the total layer was divided into 
two equal sub-layers, 500-900 mb and 100-500 mb. This took into account 
the apparent fact that balanced wind calculations in the vicinity of the 
jet stream were less reliable. 

As shown in Table 2, one of the sifnificant results was the large 
value of the flux of mean kinetic energy through the side boundaries, 
term a^, Eqn. (4) , for the upper layer, 100-500 mb. It was large enough 
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TABLE 2 

All units expressed in erg cm mb sec 
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to completely overshadow the net results both for that layer and the total 
layer, 100-900 mb. In the lower layer, 500-900 mb this dominance by- 
term a^, Eqn. (4) did not exist. In this layer the mean flux of kinetic 
energy through the upper boundary, term c^, Eqn. (4) and the eddy flux 
of kinetic energy through the side boundaries , term a^, Eqn. (4), had 
the dominant roles. For this layer, the mean and eddy terms almost 
cancelled, leaving a net -1.9 ergs an -2 mb“ ^ec - 1 . 

Results for the 200 mb jet axes were not directly comparable with 
the 300 mb axis. The balanced wind output tended to merge the two axes 
into a single axis closer to the northern axis in the first three time 
periods, and the southern axis in the last two periods. Furthermore, 
they covered only a half -wavelength. However, these results did allow 
for a check on the results presented in section IV-D of the jet stream 
at 200 mb for 16 November 1200 GMT. 

D. ANALYSIS OF AN INDIVIDUAL JET STREAM AXIS 

The results of the individual jet stream analysis are presented for 
the entrance and exit regions and are further sub-divided into a total 
layer and two sub-layers, 500-900 mb and 100-500 mb. 

Fig. 8, 9, and 10 shew the averaged vertical motion, temperature, 
and wind fields of both regions respectively. The results of the kinetic 
energy calculations for both regions are shown in Table 3. Table 3 also 
shows the result of averaging along the entire jet for use in making the 
comparison mentioned in section IV-C. 

Considering the total jet structure first, vertical motion in the 
entrance region was ascending equatorward and descending poleward of 
the jet axis (Fig. 8). The normal wind, V , was equatorward below the 

U 

jet defined at 200 mb (Fig. 10). Unfortunately, wind fields were not 
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Fig. 8 Mean vertical motion ([uil, 10 mb sec _/ ) for 
l6 November 1200 GMT 
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Fig. 9 Mean wind (m sec 1 ) and temperature (deg K) for 
l6 November 1200 GMT 
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Fig. 10 Mean perpendicular wind, v (m sec ) 
for November 1200 GMT 
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available in the balance model above the jet stream level to check if 
the transverse winds there were poleward, thus carp le ting a thermally 
direct circulation. In this region the normal winds, V^., were 
considerably stronger than those found in the exit region. 

The vertical motion in the exit region was ascending throughout, 
although much stronger poleward of the axis (Fig. 8) . The normal wind 
field, V T , did show seme poleward flow beneath the jet, but the values 
were quite small (Fig. 10) . Thus, an indication of a thermally indirect 
circulation did exist, but was not particularly well-defined. The large 
normal wind values south of the jet axis at jet stream level (Fig. 10) , 
are attributed to the location of a jet stream to the south (see Section 
II-B and Section IV-B) . 

The temperature structure (Fig. 9) , of the two regions is interesting. 
At jet stream levels, the entrance region is 3-5 K warmer than the exit 
region though the temperature patterns are the same. At 700 mb the 
temperatures for the two regions are almost identical. The large 
difference in temperatures poleward of the axis at 900 mb is attributed 
to a cold surface outbreak that occurred during this time. 

As shown in Table 3, one of the most significant results again is 
the extremely large value, for the 100-500 mb layer, of the mean and 
eddy flux of kinetic energy through the side boundaries, terms a^ and 
a^, Eqn. (4) , especially term a^. Although large values had been found 
before, none had been comparable to term a^, Eqn. (4) , as shown in 
Table 3. These results were too large to be considered reliable, and 
could not be explained by the data problems discussed in Section II-B 
and IV-B. 
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TABLE 3 



All units expressed in erg cm mb sec 
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Kinetic energy calculations using eqn. (4) 
for 16 November 1200 GET jet stream, defined 
from balanced winds at 200 mb 
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The fact that only a portion of a wavelength was analyzed, and 
that it was further sub-divided into entrance and exti regions, must be 
considered. The findings indicate results for a segment of a jet azis, 
and not the averaged results along an entire wavelength. The change of 
signs between regions. Table 3, indicates balancing between regions 
along the jet axis. Since in this case the wind acceleration is much 
stronger in the jet entrance than the deceleration in the exit zone at 
jet stream level, the regions should not be expected to completely 
balance. Notice that at lower levels, where acceleration and 
deceleration are about equal, the two regions almost balance. 

Taking the above into consideration, two conclusions result. First, 
the lower layer, 500-900 mb, is below the area where shears are large, 
and the results of the objective scheme may be considered reliable. 

Second, in the 100-500 mb layer terms g, h, and i, Eqn. (4) do not 
contain horizontal wind calculations , and term f, Eqn. (4) only uses c^. 
THus, sane comparison of these terms between layers is possible. 

In Table 3 the entire jet results are slightly different from the 
average of the two regions because of a 2° separation between entrance 
and exit regions caused when the grid was developed by the objective 
scheme. Notice that the entire jet results for the 500-900 mb layer 
are very close to those shewn in Table 2, Section IV-C, for this layer. 

Since these results are for two different situations, it is reasonable 
to assume that the results for the entrance and exit regions are reliable. 

Canparing the 500-900 mb layer jet entrance versus the jet exit 
region, it is interesting to note .the change of signs. For every term 
calculated in Eqn. (4) , except c^, the algebraic sign changed between 
regions. In addition, the entrance region shewed a net loss, while the exit 
region shewed a net gain for all terms (-16.8 vs +15.3 ergs an -2 mb- 1 sec - *) . 
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In the jet entrance region, kinetic energy was lost primarily by 
mean pressure interaction through 500 mb and the eddy flux of kinetic 
energy through the side boundaries. The exit region's kinetic energy 
gain was also primarily caused by mean pressure interaction through 
500 mb and the eddy flux of kinetic energy through the sides. But 
viewing the entire jet stream the upward transport by the mean pressure 
interaction in the entrance region is offset by downward transport in 
the exit region. 

E. INVESTIGATION OF EDDY HEAT FLUX 

This investigation was conducted on the same jet stream axis used 
in Section IV-D, and was divided into the same regions, jet entrance 
and exit zones. Fig. 11 shows the magnitudes of the terms in the heat 
budget equation, Eqn. (5), for both regions and at various levels. 

Since temperatures above 300 mb were not available from the balance 
model, only results for the region below 300 mb are shown. The vertical 
convergence of eddy heat flux and eddy adiabatic compression, term 5, 
Eqn. (5) was negligible and not included in Fig. 11. The diabatic 
temperature change, term 4, Eqn. (5) was not calculated, although sane 
estimates of it were made and are discussed later in this section. 

Mahlman (1970) suggested that the horizontal convergence of eddy 
heat flux, term 2, Eqn. (5) might counteract the expansion cooling and 
compression heating caused by the mean circulation about the jet. 

Fig. 11 shows that except for the region poleward of the jet axis at 
400 and 600 mb, eddy heat flux does not have an appreciable effect. In 
the entrance region eddy heat flux was the smallest of the four terms 
and showed a slight net transfer of heat poleward. 
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transverse vertical motion of the heat budget, eqn. 
( 5 ) for the jet stream of 16 November 1200 GKT, 
defined from balanced winds 
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In the exit region the eddy heat flux was comparable in magnitude to 
the mean horizontal transverse advection, term 1, Eqn. (5), at both 
400 and 600 mb. However, at these levels, both terms were small. 

In the jet entrance region Fig. 11 shows that the mean horizontal 
advection terms (1 and 3, Eqn. (5)) were almost balanced. The cold 
advecrtion of the mean horizontal normal advection, term 1, Eqn. (5) , 
being slightly stronger than the warm advection of the mean horizontal 
parallel advection, term 3, Eqn. (5) . In the jet exit region, the warm 
advection of the mean horizontal advection, term 3, Eqn. (5) was ■ 
approximately balanced by the expansion cooling of the mean transverse 
vertical motion, especially poleward of and at the jet axis. In the 
jet entrance this balancing of the mean horizontal advection leaves 
the mean transverse vertical motion as primarily responsible for the 
warming poleward of the jet axis in the upper levels. The mean trans- 
verse vertical motion causes cooling equatorward of the jet. In the 
jet exit region, all the terms approximately balance except at 600 and 
800 mb along the equatorward boundary, where cooling is shown. 

An interesting point shown in Fig. 11 is the reversal of roles 
between jet regions by the horizontal normal advection and the vertical 
transverse motion. Notice that in the jet entrance, the mean trans- 
verse vertical motion pattern is approximately the same as the mean 
horizontal parallel advection. In the jet exit region the mean 
horizontal normal advection pattern closely follows the mean horizontal 
parallel advection. 

The magnitudes of all terms shown in Fig. 11 decrease equatorward, 
except for the compression heating or expansion cooling due to mean 
transverse vertical motion, term 6, Eqn. (5) , in the entrance region. 
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However, as Fig. 12 shows, the local change of temperature , — — r — 
is almost constant for all three levels across the jet. Thus, it is 
implied that to maintain this structure, the magnitude of the diabatic 
tonperature change must increase equatorward. Fig. 13 shows that the 
estimated magnitude of the diabatic temperature change does increase 
equatorward, balancing the decrease equatorward of the other terms in 
Eqn. (5) . 

Confidence in the actual values of these results was limited. The 
tenperatures used were computer generated and then interpolated to the 
same levels as the balanced winds. However, fran the relative magnitude 
obtained, eddy heat flux did appear to play a role poleward of the jet 
axis. This agrees with Krishnamurti ' s (1961 b) findings for the sub- 
tropical jet that north of the jet axis there is an indication that heat 
flux is carried largely by daily eddies. The indication of this research 
is that the eddies do play a role in the heat flux, but not a dominant 
one. In addition, it appears that the diabatic tanperature change 
equatorward of the jet axis plays a significant role. 
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Fig. 12 Magnitude of the locr.1 temperature change 

eqn. (5) (in a Chr"' ) for 16 * November 1200 GMT, 
for 800 , 600 , and 400 mb 




Fig. 13 Magnitude of the estimated diabatic temperature 
change, eqn. (5) (in c Chr' / ) for l 6 November 
1200 GMT, for 800, 600, and 400 mb. 
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V. CONCLUSIONS 



Five valid conclusions are to be drawn frcm this research. First, 
the objective scheme developed did provide for a rapid and effective 
analysis of a jet stream. Wind computations in the vicinity of large 
wind shears may, however, require a finer grid than used. It is 
recommended that if wind computations are desired in large wind shear 
areas, the grid be small enough to resolve these large shears. 

Second, use of total balanced winds in the vicinity of the jet 
stream are questionable. It is recommended that observed winds be 
used above 500 mb. 

Third, there was a strong indication that the transverse circulation 
was thermally direct in the entrance zone and thermally indirect in 
the exit zone of a jet maximum. This was in agreement with the 
conclusions of Murray and Daniels (1953) . 

Fourth, kinetic energy is imported and exported frcm below the jet 
into the jet stream region primarily by the pressure interaction at 
500 mb. This expert and import along the axis is not seen if averaging 
is performed along the entire length of the axis. It is suggested that 
when studying jet stream energetics, the jet entrance and exit zones be 
included in addition to calculations along the entire length. It 
appears that seme important features in maintaining a jet maximum tend 
to be balanced between these two regions . 

Fifth, there was an indication that mean horizontal temperature 
advection and mean transverse vertical motion paly a significant role 
in the maintenance of the temperature structure poleward of the jet 
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axis. Equatorward of the jet axis there was a strong indication, from 
estimated results, that the diabatic temperature change plays a 
significant role in maintaining the thermal structure. 
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APPENDIX A 



Throughout this report frequent reference to an objective compu- 
tational schone was made. A general description of. the schone was 
given in Section II-D, but specific details were not included. The 
purpose of this Appendix is to expand on the unusual concepts and 
geometry utilized in programing. It was not designed as a step by 
step analysis of the computer program nor as a discussion of concepts 
normally utilized in programming. However, for reference, a general 
list of arguments and arrays used is listed at the end of the Appendix. 

The two concepts to be discussed are the establishment of the jet 
oriented grid system and the computing of wind components parallel and 
normal to the jet axis. In developing the concepts, for computational 
purposes, the two and a half degrees latitude and longitude between 
grid points were normalized to one. Thus, the jet axis grid points 
could be taken off the grid as if they were located on a normal 
rectangular x-y grid. Distance between grid points in the y-direction 
was 150 nautical miles. In the x-direction, distance was adjusted for 
each latitude by the formula: DX=OY *COS (PHI) , with PHI equal to the 

latitude. 

To determine the location of grid points 2.5 and 5.0 degrees lati- 
tude normal frcm the axis, the following scheme was used. First, the 
tangent at the individual jet axis grid point was determined. This 
was done by finding the distance in the x and y direction between the 
two jet axis grid points adjacent to the one used for calculation (see 
Fig. A— 1) . Using a computer library program ATAN 2 (A,B) , the angle 
T was computed (Fig. A01) . The sine and cosine of T, when added or 
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subtracted, as required, to the jet axis grid point, gives the location 
of the grid points 2 . 5 and 5 . 0 degrees latitude on either side of the 
axis. Repeating this operation at each grid point along the jet axis 
established the grid described in Section II-D. 

To calculate the wind components normal and parallel to the jet 
axis, the angle T was used in conjunction with the u and v total 
balanced wind components at the curvilinear grid points. To obtain the 
parallel component, u was multiplied by cos T and v was multiplied by 
sine T. To get the normal component, u was multiplied by sin T and v 
was multiplied by sin T (Fig. A-2) . 

To assist any research in the future that might use this program, 
the following list of arguments and arrays is included. The complete 
program may be obtained from Dr. Russell L. Elsberry, Meteorology 
Department, Naval Postgraduate School, Monterey, California 93940. 

NP = number of grid points along the jet axis. 

NC = number of grid points along the jet axis minus one. 

Used as a counter in DO-LOOPS 

NF, ND, NL, or NR = number of levels in the vertical. 

I - number of grid points along the Z or S-axis. 

K = number of grid points along the X or S-axis. 

J = number of grid points along the Y or M-axis. 

L = number of grid points along the Y or M-axis. 

TT = initial temperature field in degrees absolute. 

UU = initial zonal component of the total balanced wind in m/sec. 

W = initial meridonal component of the total balanced wind in m/sec. 

WW = initial vertical motion field in mb/sec times 10 5 . 
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GP = geopotential height fields in meters . 

PT = potential temperature. 

A C in front of the above five values means that it was the value on 
the curvilinear grid. 

PERVU = wind component normal to the jet axis. 

PARVU = wind component parallel to the jet axis. 

Any term with EDY before it designates an eddy value. All terms such 
as Al, A2, d, G, etc., refer to alphabetic designation given to various 
terms in equations discussed in the main part of this report. 
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Pig. A-l Schematic showing how angle T is calculated 
and used to define grid points 
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Fig. A-2 Schematic showing hew angle T is used 
to acquire perpendicular and parallel 
wind components 
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